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Abstract 

The cross section for the photoproduction of events containing three jets with a 
three-jet invariant mass of M^j > 50 GeV has been measured with the ZEUS detec- 
tor at HERA. The three-jet angular distributions are inconsistent with a uniform 
population of the available phase space but are well described by parton shower 
models and 0{aa^) pQCD calculations. Comparisons with the parton shower model 
indicate a strong contribution from initial state radiation as well as a sensitivity to 
the effects of colour coherence. 
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1 Introduction 



Calculations of photoproduction processes beyond leading order in perturbative QCD 
(pQCD) predict a rich variety of phenoma. Some of these can be studied in final states 
containing more than two jets. Also, the study of multijet production provides sensitive 
tests of extensions to fixed order theories such as parton shower models. The properties 
of multijet events in hadronic collisions have been the subject of earlier studies ^, ^. 
Dijet photoproduction accompanied by a third, low transverse energy cluster has been 
studied by ZEUS Q. In this paper, cross sections and angular distributions for three or 
more moderately high transverse energy jets in photoproduction are presented for the 
first time. 

Apart from the azimuthal orientation, a system of two massless jets can be completely 
specified in its centre-of-mass (CM) frame by the two-jet invariant mass, M2J and cos-d*, 
where is the angle between the jet axis and the beam-line. The distribution in cos'f?* 
for photoproduction of dijets is forward-backward peaked with sensitivity to the spin of 
the exchanged fermion or boson 

A set of observables describing events with an arbitrary number of jets has been proposed 
which spans the multijet parameter space, facilitates the interpretation of the data within 
pQCD and reduces to M23 and cos??* for the dijet case Q. For three massless jets there 
are five parameters which are defined in terms of the energies, Ei, and momentum 
three-vectors, pi, of the jets in the three-jet CM frame and pB, the beam direction]^ 
The jets are numbered, 3, 4 and 5 in order of decreasing energy as illustrated in the 
schematic drawing. Fig. |l|. The parameters are: the three-jet invariant mass, M^j; the 
energy-sharing quantities X3 and X4, 

the cosine of the scattering angle of the highest energy jet with respect to the beam, 

cosi?3 = — — — ; (2) 

\pb\\P3\ 

and the angle between the plane containing the highest energy jet and the beam 
and the plane containing the three jets. The latter is defined by 

, im X Pb) ■ (P4 X ps) 
cos 1^3 = — — — — . (3) 

|P3 X PBim X P5\ 

The definition of the angles 'd^ and -03 is illustrated in Fig. |^. Since 'ds involves only the 
highest energy jet, the distribution of cos 1)3 in three-jet processes may be expected to 
follow closely the distribution of cos tD* in dijet events. The ip^ angle, on the other hand, 
reflects the orientation of the lowest energy jet. In the case where this jet arises from 
initial-state radiation, the coherence property of QCD will tend to orient the third jet 
close to the incoming proton or photon direction. The two planes shown in Fig. 1 will 
therefore tend to coincide leading to a ip^ distribution which peaks toward and vr. 

^We take the nominal beam direction as ps ~ z. In the ZEUS coordinate system the z-axis is 
defined to be in the proton beam direction. Polar angles, §, are measured with respect to the z-axis 
and pseudorapidity is defined as 77 = — ln(tan -1). 
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This paper presents the three-jet inclusive cross section in photoproduction and the 
distribution of the three-jet events with respect to M3J, X3, X4, cos-i^s and tp^. This 
work was performed with the ZEUS detector using 16 pb~^ of data dehvered by HERA 
in 1995 and 1996. 



2 Experimental Conditions 

In this period HERA operated with protons of energy Ep = 820 GeV and positrons of 
energy = 27.5 GeV. The ZEUS detector is described in detail in @, |8|. The main 
components used in the present analysis are the central tracking system positioned 
in a 1.43 T solenoidal magnetic field and the uranium-scintillator sampling calorimeter 
(CAL). The tracking system was used to establish an interaction vertex. Energy deposits 
in the CAL were used in the jet finding and to measure jet energies. The CAL is 
hermetic and consists of 5918 cells each read out by two photomultiplier tubes. Under 
test beam conditions the CAL has energy resolutions of 18%/y/E (GeV) for electrons 
and 35%/ a/ (GeV) for hadrons. Jet energies were corrected for the energy lost in 
inactive material in front of the CAL which is typically about one radiation length 



(see Section 3.4). The effects of uranium noise were minimized by discarding cells 
in the electromagnetic or hadronic sections if they had energy deposits of less than 
60 MeV or 110 MeV, respectively. The luminosity was measured from the rate of the 
bremsstrahlung process e~^p — > e'^pj. A three-level trigger was used to select events 
online . 



3 Analysis 

3.1 Offline Cleaning Cuts 

To reject residual beam-gas and cosmic ray backgrounds, tighter cuts using the final 
2-vertex position, other tracking information and timing information are applied of- 
fline. Neutral current deep inelastic scattering (DIS) events with an identified scattered 
positron candidate in the CAL are removed from the sample as described in detail else- 
where p. Charged current DIS events are rejected by a cut on the missing transverse 
momentum measured in the CAL. Finally, a restriction is made on the range of y, the 
fraction of the positron's energy carried by the incoming photon. The requirement, 
0.15 < yjB < 0.65 is made where yjB is an estimator of y which is determined from the 
CAL energy deposits according to the Jacquet-Blondel method |jl^. This requirement 
corresponds to accepting events in the range 0.2 < y < 0.8. These cuts restrict photon 
virtualities to less than about 1 GeV^ with a median of around 10~^ GeV^. 



3.2 Jet Finding 



Jets are found using the KTCLUS |11] finder in the inclusive mode [12|. This is a 
clustering algorithm which combines objects with small relative transverse energy into 
jets. It is invariant under Lorentz boosts along the beam axis and is ideal for the study of 
multijet processes since it suffers from no ambiguities due to overlapping jets. Once the 
jets are determined, their transverse energy, pseudorapidity and azimuth are calculated 

4'* = Et,, r?j-t = (1/4*) Ei ET.Vi and 



according to the Snowmass convention [O 
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(pi^^ = J2i Exi^i, where the sum runs over all objects assigned to the jet. The 

energies and three- vectors of the jets are then determined from the Ej^^, r^^* and ip^^^. 
The objects input to the jet algorithm may be hadrons in a simulated hadronic final 
state, the final state partons of a pQCD calculation, or energy deposits in the detector. 
In the following, a jet quantity constructed from CAL cells with no energy correction 
has the superscript "CAL" while a jet quantity constructed from CAL cells and then 
subjected to a correction for energy loss in inactive material has the superscript "COR". 
There is no additional superscript for quantities referring to jets of final state partons 
or hadrons. 



3.3 Monte Carlo Event Simulation 

The response of the detector to jets and the acceptance and smearing of the measured 
distributions are determined using samples of events generated from Monte Carlo (MC) 
simulations. We have used the programs PYTHIA 5.7 Q and HERWIG 5.9 
which implement the leading order matrix elements followed by parton showers. In these 
simulations multijet events can originate through this parton shower mechanism. Colour 
coherence in the parton shower is treated differently in the two models. In PYTHIA, 
parton showers are evolved in the squared mass of the branching parton with colour 
coherence effects implemented as a restriction on the opening angle of the radiation. In 
contrast, in HERWIG a parton shower evolution variable is chosen which automatically 
limits the branching to an angular ordered region. For both models, leading order direct 
and resolved processes are generated separately and combined according to the ratio of 
their generated cross sections. For the uncorrected distributions presented in this section 
the minimum transverse momentum of the partonic hard scatter (p™™) was set to 4 GeV. 
In Section ^, corrected cross sections are presented and compared with the predictions 
from HERWIG and PYTHIA with p™™ = 8 GeV (our conclusions are insensitive to this 
parameter). The photon parton densities used were GRV LO and the proton parton 



densities were CTEQ4 LO [17|. In the HERWIG simulation of the resolved processes. 



multiparton interactions have been included (these are not important in this kinematic 
regime, as discussed in Section Q). 

The quality of these simulations is illustrated in Fig. ^ which shows the transverse 
energy fiow around the jet axes. In this comparison the events have two jets with 
^jet CAL ^ g GeV and a third jet with E';^* > 4 GeV and jet pseudorapidities 
jj^et CAL I ^ 2.4. The additional requirements on the CM quantities, M£^^ > 42 GeV, 
I cost^al*^^^ < 0.8 and X^^^ < 0.95, have also been applied. These conditions represent 
those on the selected events (described in Section |3.5| ) to a good approximation. The 
jets have a narrow core with little transverse energy in the pedestal, except for the 
lowest E}^^ jet where significant contribution to the "pedestal" from the other two 
jets in the event would be expected. In the A(f profiles of the two highest i?;^* '^^^ 
jets, the peaks near ibvr indicate that these are roughly back-to-back. The PYTHIA 
and HERWIG event samples were passed through a detailed simulation of the ZEUS 
detector and the same selection criteria as for the data were applied. The MC samples 
provide a reasonable description of the energy flow in these three-jet events. HERWIG 
generates somewhat too little transverse energy in regions far from the jet core in ip for 
the lowest transverse energy jet, while PYTHIA slightly overestimates the transverse 
energy in the core of this jet. These models are also able to reproduce satisfactorily the 
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yjB distribution for the three-jet events, the lab- frame -Ej'* rf^^ '^^^ distributions 

as weh as the transverse and longitudinal components of the boost from the lab-frame 
to the CM- frame (not shown). 

3.4 Jet Energy Corrections 

Jets of hadrons lose about 15% of their transverse energy when passing through inactive 
material before impinging on the CAL. This energy loss has been corrected using the 
MC samples ^]. The KTCLUS algorithm was applied to the hadronic final state 
and from comparison of these hadron jets with the CAL jets obtained after the detector 
simulation, correction factors were determined as a function of E^r^^ and r/j'^* CAL^ 
After applying these corrections the average shift in Ei^^ within the MC simulation is 
less than 2% and the Ej^^ resolution is 14%. This may be compared with the global 
jet energy scale uncertainty of ±5% |^]. The correction also reduces the shift in the 
reconstruction of M3J from 16% to less than 1%. After these corrections for jet energy 
loss the average resolutions are 8% in M3J, 0.03 units in X3, 0.05 units in X4, 0.03 units 
in COSIES and 0.1 radians in ^3 and the distributions are well centred on their expected 
values. 

3.5 Event Selection 

After the jet energy correction the events are required to have at least two jets with 
COR ^ g GeV, a third jet with E^^^ > 5 GeV and jet pseudorapidities in the 
range Irf'^^ COR| ^ 2.4. The requirement of high transverse energy for the jets ensures 
that the process should be calculable within pQCD. However, it introduces a bias in the 
angular distributions by excluding jets that are produced close to the beam-line. We 
make the additional requirements M£°^ > 50 GeV, | cos < 0.8 and X^'^^ < 0.95 

to minimize such a bias. After these cuts the mean transverse energy of the highest, 
second-highest and third-highest transverse energy jet is about 20 GeV, 15 GeV and 
10 GeV. From 16 pb^^ of data, 2821 events are selected. Around 15% of these have a 
fourth jet with Eij?* > 5 GeV, in agreement with the prediction of the parton shower 
models. Backgrounds from beam gas and cosmic ray events, determined from unpaired 
bunch crossings, are negligible. The DIS contamination, determined using Monte Carlo 
techniques, is around 1% and neglected. 

3.6 Acceptance Correction 

The MC samples have been used to correct the data for the inefficiencies of the trig- 
ger and the offline selection cuts and for migrations caused by detector effects. The 
correction factors are calculated as the ratio Ntme/^rec in each measured bin where 
Ntrue is the number of events generated in the bin and A'rec is the number of events 
reconstructed in the bin after detector smearing and all experimental cuts. The final 
bin-by-bin correction factors lie between about 0.7 and 1.3, the dominant effect arising 
from migrations across the M3J threshold. The cross sections were determined using the 
corrections obtained with PYTHIA. 
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3.7 Systematic Uncertainties 

A detailed study of the sources contributing to the systematic uncertainties of the mea- 
surements was performed ||T^. Only the significant sources are listed here. 

• The acceptance correction was performed using HERWIG instead of PYTHIA. 
The uncertainties associated with the model are typically around 20% and this 
forms the dominant uncertainty on the area-normalized distributions. 

• The absolute energy scale of the detector response to jets with Ej^^ > 5 GeV is 
known to ±5% This leads to an uncertainty of 15 to 20% on the cross section. 
This is the dominant systematic uncertainty on the normalization of the cross 
section but as this uncertainty is highly correlated between bins it has a negligible 
affect on the area- normalized distributions. 

• The results were recalculated allowing for fluctuations from outside the selected 
kinematic region by relaxing each of the cut parameters by 1 a of the resolution. 
This effect is typically 5%. 

The systematic uncertainties have been added in quadrature to the statistical errors and 
this is shown as the outer error bars in the figures, with the exception of the absolute 
jet energy scale uncertainty which is shown as a shaded band for the cross sections. An 
overall normalization uncertainty of 1.5% from the luminosity determination has not 
been included. 



4 Results and Discussion 

The three-jet inclusive cross section is presented for events having at least two jets with 
> 6 GeV and a third jet with > 5 GeV where the jets satisfy \r^''^\ < 2.4. 
This cross section refers to jets in the hadronic final state. To minimize the effects of 
these jet cuts on the distributions of physical interest, the requirements M3J > 50 GeV, 
|cosi?3| < 0.8 and X3 < 0.95 have been imposed. The cross section is presented for 
photon-proton CM energies W-yp in the range 134 GeV< W-yp < 269 GeV and the 
negative square of the invariant mass of the incoming photon extending to 1 GeV^ . The 
cross section is cr = 162 ± 4(stat.)t6^(sys.)lt25(6nergy scale) pb. 

A study using the PYTHIA MC indicates that hadronization effects are small (~ 5%), 
and flat in the distributions presented here [19|. The measurements are directly con- 



fronted with 0{aa1) pQCD calculations from two groups of authors |20, 21, 22]. The 
CTEQ4 LO proton parton densities and the GRV LO J^] photon parton densities 
have been used in these calculations. The renormalization and factorization scales, fi, 
have been chosen to equal Ejf^^, where Ej^^^ is the largest of the Ej^^ values of the 

(5) 

three jets. was calculated at one loop with A^ = 181 MeV. As the calculations are 
leading order for three-jet production the normalization uncertainty due to the choice 
of /i is expected to be large. An uncertainty of a factor of two in the cross section for 
variation of fi between E^^^/2 and 2Ej^^^ has been quoted ]23]. 

The three-jet invariant mass distribution is shown in Fig. |. The cross section falls 
approximately exponentially from the threshold value at 50 GeV to the highest measured 
value, around 150 GeV. The data are compared with the two 0{aa'^) pQCD calculations. 



5 



These are in good agreement with the data, even though the calculations are leading 
order for this process. The M3J distributions predicted by the parton shower models 
PYTHIA and HERWIG are also in agreement with the data in shape although the 
predicted cross sections are too low by 30-40%. 

The distributions of the fraction of the available energy taken by the highest and second- 
highest energy jets are shown in Figs. ^ (a) and (b), respectively. Here the prediction for 
three jets uniformly distributed in the available phase space (i.e. with a constant matrix 
element) is also shown as the dotted curve. The parton shower models give a reasonable 
description of these energy sharing quantities. The pQCD calculations (overlapping) 
are in excellent agreement with these distributions. However, the similarity between the 
measured distributions and the three-body phase space prediction indicates that these 
distributions have little sensitivity to the pQCD matrix elements. 

In Figs. ^ (c) and (d) the cosies and tp^ distributions are shown. These angular dis- 
tributions are dramatically different from the distributions obtained from phase space, 
demonstrating that these quantities are sensitive to the pQCD matrix elements. The 
cos??3 distribution has forward and backward peaks, as expected. The O(aa^) pQCD as 
well as the parton shower calculations, which take into account the dependence of this 
distribution on the spin of the exchanged quark or gluon, are in good agreement with the 
data. The jet algorithm and minimum i^j?* requirements deplete the data near 1/^3 ~ 
and vr as indicated by the shape of the phase space curve in Fig. |(d). With this taken 
into account, the data indicate a strong tendency for the three-jet plane to lie near 
the plane containing the beam and the highest energy jet. This effect is reproduced 
in the 0(q:q;^) matrix element calculations. It is interesting that the parton shower 
Monte Carlo programs PYTHIA and HERWIG are also able to provide a reasonable 
representation of the shape of the ips distribution. 

Including a simulation of multiparton interactions in the parton shower programs has 
been found to improve significantly the description of low E}^^ photoproduction Q . In 
the present study the sensitivity to multiparton interactions has been investigated using 
both PYTHIA and HERWIG |19|. In neither case do secondary parton interactions 
cause a significant difference in the three-jet cross section in this kinematic regime, or in 
the angular distributions generated. It appears therefore that the third jet arises here 
from the parton shower and is not due to a second hard scatter. 

Within the parton shower prescription, it is possible to separate the contributions of 
initial and final state parton showers. In Fig. ^(a) the three-jet cross section as a function 
of '03 is shown and compared with the predictions of PYTHIA. The MC events have 
been separated into three samples; initial-state radiation only, final-state radiation only, 
and default PYTHIA which includes the interference of these two. The area-normalized 
distributions of ip^ are compared with these models in Fig. |5|(b). Both the normalization 
and shape of these distributions indicate that the observed three-jet production occurs 
predominantly through initial state radiation with the final state radiation making a 
small contribution. 

The QCD phenomenon of colour coherence is implemented in the PYTHIA parton 
shower model by prohibiting radiation into certain angular regions which are determined 
by the colour flow of the primary scatter. It is possible within this model to switch 
QCD colour coherence on and off. Figs, ^(c) and (d) again show the cross section and 
the area-normalized distribution with respect to V's compared with the HERWIG and 
PYTHIA predictions. The data lie above these predictions as previously mentioned. 
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however this discrepancy is not regarded as significant in view of the limited order 
of the calculation. The predictions do reproduce reasonably well the shape of the ips 
distribution. This is not the case if the simulation is done with colour coherence switched 
off. The incoherent PYTHIA prediction is relatively flat in ip^. Coherence reduces the 
phase space available for large angle emissions as indicated by the drop in cross section 
around V's ~ 7r/2 for default PYTHIA. Colour coherence in the parton shower is needed 
to describe the shape of this distribution. QCD colour coherence seems to be a stronger 
effect in HERWIG than in PYTHIA however the present data are not precise enough 
to discriminate between these two simulations. 

5 Summary 

The inclusive cross section for the photoproduction of three jets has been measured 
by the ZEUS collaboration at HERA. 0{aa1) pQCD calculations are able to describe 
the cross section da/dM^j, while parton shower models underestimate the cross section 
but are consistent in shape with the M3J distribution. The angular distributions of 
the three jets are inconsistent with a uniform population of the available phase space 
but are well described by both fixed-order pQCD calculations and parton shower Monte 
Carlo models. Simulation of multiparton interactions does not help to describe the 
data in this kinematic regime. Within the parton shower model the three-jet events 
are found to occur predominantly due to initial state radiation, and the fundamental 
QCD phenomenon of colour coherence is seen to have an important effect on the angular 
distribution of the third jet. 
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THREE-BODY REST FRAME 




Figure 1: Illustration of the angles ■i?3 and -03 for a particular three-jet configuration. 
The beam direction is indicated hy pB- 
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Figure 2: Uncorrected transverse energy flow with respect to the jet axis in the labora- 
tory frame for three-jet events in order of -B^* CAL j^£^ uncorrected energy 
flow with respect to cp is shown for cells within one unit of r/ of the jet axis while on the 
right the profile with respect to rj is shown for cells within one radian of ip of the jet 
axis. The data arc shown as black dots while the PYTHIA and HERWIG predictions 
are shown by the solid and dashed histograms, respectively. 
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Figure 3: The measured three-jet cross section with respect to the three-jet invariant 
mass, da/dM^j, is shown by the black dots where the inner error bar shows the statis- 
tical error and the outer error bar is the sum in quadrature of the statistical error and 
the systematic uncertainty. The jet energy-scale uncertainty, which is highly correlated 
between bins, is shown separately as the shaded band. O(aa^) pQCD calculations by 
Harris & Owens and Klasen are shown by the thick solid and dot-dashed lines, respec- 
tively. The thin solid and dashed histograms show the predictions from two different 
parton shower models, PYTHIA and HERWIG. 
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Figure 4: The distributions of the energy sharing quantities, and X^, are shown by 
the black dots in (a) and (b) , respectively, and the distributions of the cos t?3 and V's are 
shown in (c) and (d) . Inner error bars show the statistical error and the outer error bars 
show the quadratic sum of this with the systematic uncertainty. The fixed-order pQCD 
predictions are shown by the thick solid and dot-dashed lines and the parton shower 
model predictions are shown by the thin solid and dashed histograms. The phase space 
distribution of three jets is indicated by the dotted line. 
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Figure 5: The measured cross section da/dip^ is shown in (a) and (c) and the area- 
normahzed distribution of tps is shown in (b) and (d). The error bars are as described 
previously with the correlated systematic uncertainty due to the jet energy-scale shown 
as the shaded band in (a) and (c). The solid histogram shows the default PYTHIA 
prediction. In (a) and (b) the dashed and dot-dashed histograms show the predictions 
from PYTHIA with final state radiation switched off and with initial state radiation 
switched off. In (c) and (d) the dashed and dot-dashed histograms show the predictions 
of HERWIG and of PYTHIA with colour coherence switched off. 
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